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Abstract 

o 

04 ' Physarum Polycephalum is a slime mold that is apparently able to solve shortest path 

problems. A mathematical model has been proposed by biologists to describe the feedback 
mechanism used by the slime mold to adapt its tubular channels while foraging two food 
sources sq and si. We prove that, under this model, the mass of the mold will eventually 
■ converge to the shortest sq-si path of the network that the mold lies on, independently 

of the structure of the network or of the initial mass distribution. 

This matches the experimental observations by the biologists and can be seen as an 
I example of a "natural algorithm", that is, an algorithm developed by evolution over 

. millions of years. 
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Figure 1: The experiment in |NYT00j (reprinted from there): (a) shows the maze uniformly 
covered by Physarum; the yellow color indicates the presence of Physarum. Food (oatmeal) 
is provided at the locations labelled AG. After a while, the mold retracts to the shortest path 
connecting the food sources as shown in (b) and (c). (d) shows the underlying abstract graph. 
The video [YoulOj shows the experiment. 

1 Introduction 

Physarum Polycephalum is a slime mold that is apparently able to solve shortest path prob- 
lems. Nakagaki, Yamada, and Toth |NYTOO| report on the following experiment, see Figured) 
They built a maze, covered it with pieces of Physarum (the slime can be cut into pieces that 
will reunite if brought into vicinity), and then fed the slime with oatmeal at two locations. 
After a few hours, the slime retracted to a path that follows the shortest path connecting the 
food sources in the maze. The authors report that they repeated the experiment with dif- 
ferent mazes; in all experiments, Physarum retracted to the shortest path. There are several 
videos available on the web that show the mold in action [You 10] . 

Paper |TKN07j proposes a mathematical model for the behavior of the slime and argues 
extensively that the model is adequate. We will not repeat the discussion here but only define 
the model. Physarum is modeled as an electrical network with time varying resistors. We 
have a simple undirected graph G = {N,E) with distinguished nodes sq and si, which model 
the food sources. Each edge e & E has a positive length Lg and a positive diameter De{t); 
Lf. is fixed, but D(,{t) is a function of time. The resistance Re{t) of e is Re{t) = L(./De(t). 
We force a current of value 1 from sq to si. Let Qe{t) be the resulting current over edge 
e = {u,v), where {u,v) is an arbitrary orientation of e. The diameter of any edge e evolves 
according to the equation 



flow through any edge is equal to its diameter. In non-equilibrium, the diameter grows or 
shrinks if the absolute value of the flow is larger or smaller than the diameter, respectively. In 
the sequel, we will mostly drop the argument t as is customary in the treatment of dynamical 
systems. 

The model is readily turned into a computer simulation. In an electrical network, every 
vertex v has a potential is a function of time. We may fix ps-^ to zero. For an edge 

e = {u,v), the flow across e is given by {pu —pv)/Re- We have flow conservation in every 
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vertex except for sq and si; we inject one unit at sq and remove one unit at si. Thus, 



where 6{v) is the set of neighbors of v and 6^0 = 1; ^si = ~1) and 6^ = otherwise. The 
node potentials can be computed by solving a linear system (either directly or iteratively) . 
Tero et al. |TKN07j were the first to perform simulations of the model. They report that 
the network always converges to the shortest sq-si path, i.e., the diameters of the edges on 
the shortest path converge to one, and the diameters on the edges outside the shortest path 
converge to zero. This holds true for any initial condition and assumes the uniqueness of the 
shortest path. 

Miyaji and Ohnishi |MQ07|. IMOOSj initiated the analytical investigation of the model. 
They argued convergence against the shortest path if G is a planar graph and sq and si lie 
on the same face in some embedding of G. 

Our main result is a convergence proof for all graphs. For a network G = (V, E, so,si,L), 
where {Le)e£E is a positive length function on the edges of G, we use Go = {V, Eq) to denote 
the subgraph of all shortest source-sink paths, L* to denote the length of a shortest source- 
sink path, and £* to denote the set of all source-sink flows of value one in Gq. If we define 
the cost of flow Q as ^g-LgQej then £* is the set of minimum cost source-sink flows of value 
one. If the shortest source-sink path is unique, £* is a singleton. The dynamics are attracted 
by a set A C M.^ if the distance (measured in any Lp-norm) between D{t) and A converges 
to zero. 

Theorem [Theorem [2] in Section [6] Let G = {V, E, sq, si, L) be an undirected network 
with positive length function {Le)e^E- Let -De(O) > 6e the diameter of edge e at time zero. 
The dynamics ([7P are attracted to 8* . If the shortest source-sink path is unique, the dynamics 
converge to the flow of value one along the shortest source-sink path. 

We conjecture that the dynamics converge to an element of S* but only show attraction 
to £* . A key part of our proof is to show that the function 

V = . ^ „ + (^{«o} - 1)' (3) 

decreases along all trajectories that start in a non-equilibrium configuration. Here, C is the 
set of all sq-si cuts, i.e., the set of all 5 C iV with sq G S and si ^ S; Cs = J2e€5{S) 
is the capacity of the cut S when the capacity of edge e is set to Df,; and min^gc Cs (also 
abbreviated by C) is the capacity of the minimum cut. The first term in the definition of V 
is the normalized hardware cost; for any edge, the product of its length and its diameter may 
be interpreted as the hardware cost of the edge; the normalization is by the capacity of the 
minimum cut. We will show that the first term decreases except when the maximum flow 
F in the network with capacities Df. is unique, and moreover, \Qe\ = \F(.\/C for all e. The 
second term decreases as long as the capacity of the cut defined by sq is different from 1. We 
show that the capacity of the minimum cut converges to one and that the derivative of V is 
upper bounded by — X]e(-^min/4)(-De/C' — IQel)^, where Lmm is the minimum length of any 
edge. Since V is non-negative, this will allow us to conclude that \D(, — \Qe\ \ converges to zero 
for all e. In the next step, we show that the potential difference A = p^,, ~Psi between source 
and sink converges to the length L* of a shortest-source sink path. We use this to conclude 
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Figure 2: Part (a) illustrates the path decomposition. All edges are assumed to have length 
1; Po = (ei), Pi = (62,63,64), P2 = (65,66), Pl = 1, pt, =0,pl = 1/3, pi = 2/3, p*^ = 1/2, 
/(Pi) = 1/3, and /(P2) = 1/6. 

Part (b) shows the Wheatstone graph. The direction of the flow on edge {u, v} may change 
over time; the flow on all other edges is always from left to right. 

that De and Qe converge to zero for any edge e ^ Eq. Finally, we show that the dynamics 
are attracted by 8* . 

We found the function V by analytical investigation of a network of parallel links (see 
Section H]), extensive computer simulations, and guessing. Functions decreasing along all 
trajectories are called Lyapunov functions in dynamical systems theory |HS74| . The fact that 
the right-hand side of system ([1]) is not continuously differentiable and that the function V 
is not differentiable everywhere introduces some technical difficulties. 

The direction of the flow across an edge depends on the initial conditions and time. We 
do not know whether flow directions can change infinitely often or whether they become 
ultimately fixed. Under the assumption that fiow directions stabilize, we can characterize 
the (late stages of the) convergence process. An edge e = {n, v} becomes horizontal if 
limt^oo \Pu — Pv\ =0, and it becomes directed from n to f (directed from v to u) if pu > Pv 
for all large t {p^ > pu for all large t). An edge stabilizes if it either becomes horizontal or 
directed, and a network stabilizes if all its edges stabilize. If a network stabilizes, we partition 
its edges into a set Eh of horizontal edges and a set ^ of directed edges. If {u,v} becomes 
directed from u to then (n, v) € 

We introduce the notion of a decay rate. Let r < 0. A quantity D{t) decays with rate at 
least r if for every e > there is a constant A such that InP'(t) < A + {r + e)t for all t. A 
quantity D{t) decays with rate at most r if for every e > there is a constant a such that 
\aD{t) > a + (r — £)t for all t. A quantity D[t) decays with rate r if it decays with rate at 
least and at most r. 

Lemma [Lemma 1201 in Section [7j For e € Eh, decays with rate —1 and \Qe\ decays 
with rate at least —1. 

We define a decomposition of G into paths Pq to Pk, an orientation of these paths, a slope 
/(P.) for each P., a vertex labeUing p*, and an edge labehmg r. Po is 43 shortest sq-si path in 
G, /(Pq) = 1, Te = /(Po) — 1 for all e G Pq, and = dist(t>, si) for all v € Pq, where dist(t>, si) 

^We assume that Po is unique. 
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is the shortest path distance from v to si. For I < i < k, we hav^ Pj = argmaxpgp fiP), 
where V is the set of all paths P in G with the following properties: (1) the startpoint a 
and the endpoint 6 of P lie on Pq U . . . U Pi_i, > pi, and /(P) = (p* - pl)/L{P); (2) no 
interior vertex of P lies on Pq U . . . U Pi_i; and (3) no edge of P belongs to Pq U . . . U Pj_i. 
If Pa > Pb^ direct Pj from a to 6. If p* = pi, we leave the edges in Pj undirected. We set 
re = f{Pi) — 1 for all edges of Pi, and p% = pI + f{Pi) distp^(t', b) for every interior vertex v 
of Pj. Figure [2{a) illustrates the path decomposition. 

Lemma [Lemma 1211 in Section \7\ There is an iq < k such that 

/(Po) > /(Pi) > . . . > f{Pi,) > = /(P^o+i) = . . . = /(Pfc). 



Theorem [Theorem [3] in Section \7\ If a network stabilizes, — Uj<jgP(Pj), the orien- 
tation of any edge e G ^ agrees with the orientation induced by the path decomposition, and 
Eh = Ui>i(,P(Pj). The potential of each node v converges to p*. The diameter of each edge 
e G P \ Pq decays with rate r^ . 

We cannot prove that flow directions stabilize in general. For series-parallel graphs, flow 
directions trivially stabilize. The Wheatstone graph, shown in Figure [2]Jb), is the simplest 
graph, in which flow directions may change over time. 

Theorem [Theorem [5] in Section [5] The Wheatstone graph stabilizes. 

The uncapacitated transportation problem generalizes the shortest path problem. With 
each vertex v, a supply/demand by is associated. It is assumed that Ylv^v = 0. Nodes with 
positive by are called supply nodes, and nodes with are called demand nodes. 

In the shortest path problem, exactly two vertices have non-zero supply /demand. A feasible 
solution to the transportation problem is a flow / satisfying the mass balance constraints, i.e., 
for every vertex v, by is equal to the net flow out of v. The cost of a solution is Lef^. The 
Physarum solver for the transportation problem is as follows: At any fixed time, the potentials 
are defined by ([2]) and the currents {Qe)eeE are derived from the potentials by Ohm's law. 
The dynamics evolve according to ([T]). The equilibria, i.e., \Qe\ = for all e, are precisely 
the flows with the following equal-length property. Orient the edges in the direction of Q and 
drop the edges of flow zero. In the resulting graph, any two distinct directed paths with the 
same source and sink have the same length. Let £ be the set of equilibria. 

Theorem [Theorem [8] in Section [9] The dynamics ^^ are attracted to the set of equilibria 
8. If any two equilibria have distinct cost, the dynamics converge to an optimum solution of 
the transportation problem. 

The convergence statement for the transportation problem is weaker than the correspond- 
ing statement for the shortest path problem in two respects. There, we show attraction to 
the set of equilibria of minimum cost (now only to the set of equilibria) and convergence to 
the optimum solution if the optimum solution is unique (now only if no two equilibria have 
the same cost). 

This paper is organized as follows: In Section [2l we discuss related work, and in Section [Sj 
we put our results into the context of natural algorithms and state open problems. The 

■^We assume that Pi is unique except if f{Pi) — 0. 
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technical part of the paper starts in Section [H We first treat a network of parallel links; this 
situation is simple enough to allow an analytical treatment. In Section [SJ we review basic 
facts about electrical networks and prove some simple facts about the dynamics of Physarum. 
In Section we prove our main result, the convergence for general graphs. In Section [TJ 
we prove exponential convergence under the assumption that flow directions stabilize, and 
in Section [8l we show that the Wheatstone network stabilizes. Finally, in Section [9l we 
generalize the convergence proof to the transportation problem. 



2 Related Work 

Miyaji and Ohnishi [MO071 IMQOSj initiated the analytical investigation of the model. They 
argued convergence against the shortest path if G is a planar graph and sq and si lie on the 
same face in some embedding of G. Ito et al. |IJNTll| study the dynamics ([T|) in a directed 
graph G = {V, E); they do not claim that the model is justified on biological grounds. Each 
directed edge e has a diameter D(.- The node potentials are again defined by the equations 

The summation on the right-hand side is over all neighbors u of v; edge directions do not 
matter in this equation. If there is an edge from n to f and an edge from v to u, u occurs 
twice in the summation, once for each edge. The dynamics for the diameter of the directed 
edge uv are then Duv = Quv — D^v, where Quv = Duvipu — Pv)/Luv The dynamics of this 
model are very different from the dynamics of the model studied in our paper. For example, 
assume that there is an edge vu, no edge uv, and pu > Pv always. Then Q^u < always and 
hence D^u will vanish at least with rate —1. The model is simpler to analyze than our model. 
Ito et al. prove that the directed model is able to solve transportation problems and that the 
Dg's converge exponentially to their limit values. 



3 Discussion and Open Problems 

Physarum may be seen as an example of a natural computer, i.e., a computer developed by 
evolution over millions of years. It can apparently do more than compute shortest paths and 



solve transportation problems. In TTS^IO , the computational capabilities of Physarum are 
applied to network design, and it is shown in lab and computer experiments that Physarum 
can compute approximate Steiner trees. No theoretical analysis is available. The book [AdalO] 
and the tutorial [NTK+09] contain many illustrative examples of the computational power of 
this slime mold. 

Chazelle | Cha09j advocates the study of natural algorithms; i.e., "algorithms developed 
by evolution over millions of years", using computer science techniques. Traditionally, the 
analysis of such algorithms belonged to the domain of biology, systems theory, and physics. 
Computer science brings new tools. For example, in our analysis, we crucially use the max- 
flow min-cut theorem. Natural algorithms can also give inspiration for the development of new 
combinatorial algorithms. A good example is |CKM"'"ll] . where electrical flows are essential 
for an approximation algorithm for undirected network flow. 

We have only started the theoretical investigation of Physarum computation, and so many 
interesting questions are open. We prove convergence for the dynamics = f{\Qe\) — D^, 
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where / is the identity function. The biological literature also suggests the use of /(x) = 
x'^ /{I + x^) for some parameter 7. Can one prove convergence for other functions /? We 
prove that flow directions stabilize in the Wheatstone graph. Do they stabilize in general? 
We prove, but only for stabilizing networks, that the diameters of edges not on the shortest 
path converge to zero exponentially for large t. What can be said about the initial stages of 
the process? The Physarum computation is fully distributed; node potentials depend only 
on the potentials of the neighbors, currents are determined by potential differences of edge 
endpoints, and the update rule for edge diameters is local. Can the Physarum computation 
be used as the basis for an efficient distributed shortest path algorithm? What other problems 
can be provably solved with Physarum computations? 

4 Parallel Links 

We discovered the Lyapunov function used in the proof of our main theorem through exper- 
imentation. The experimentation was guided by the analysis of a network of parallel links. 
In such a network, there are vertices sq and si connected with m edges of lengths Li < L2 < 
. . . < L^. Let Di be the diameter of the i-th link, and let D = Di. Let A = — 
be the potential difference between source and sink. Then, Qi = A/Ri = DiA/Li. Since 
Qi = 1, we have A = 1/ A/i^i- 

Lemma 1 The equilibrium points are precisely the single links. 

Proof: In an equilibrium point, Qi = Di for all i. Since Qi = DiA/Li, this implies A = Lj 
whenever Qi / 0. Thus, in an equilibrium there is exactly one i with Qi 7^ 0. Then, Qi = 1. 



Lemma 2 Let D = Di. Then, D converges to 1. 

Proof: have D = Di = Qi — Di = 1 — D. The claim follows by directly solving 
the differential equation: D{t) = 1 + {D{0) - l)exp(-t). I 

For networks of parallel links, there are many Lyapunov functions. 

Lemma 3 Let Xi = Di/D, and let L be such that 1/L = Xj/Lj. The quantities 

'^Di/D, ^XjLj, L, ^QjLj, A^ALj, and ^(Lj In A - -^1 In -Di) 

j>2 i i i i>2 

decrease along all trajectories, starting in non- equilibrium points. 

Proof: Clearly, J2j = 1 ^^d A = L/D. The derivative Xi of Xi computes as: 

. _ DjD - DjP _ {DiA/Li-Di)D-Di{l-D) _ f _ l\ _ ]_ f ^ _ ,\ 
" £»2 - D2 - [l.d dJ'''~ D\Li J 

We have L > Li iff ^j>2 > 0- Thus, the derivative of xi is zero if xi = 1 and positive if 
xi < 1. Thus, X]j>2 ^2 decreases along all trajectories, starting in non-equilibrium points. 
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Let V = YliiXiLi. Then, 



So, it suffices to show ^-LiXi > L = l/Yli^i/Li, or equivalently, iJ2iLiXi){^iXi/Li) > 1. 
This is an immediate consequence of the Cauchy-Schwarz inequahty. Namely, 

Now, let y = 1/L = Xj/Lj. We show that V is increasing. We have 

LXi 1 X. 



Tji D y Li J Li D y Li Li Li 



Let Zi = Lxi/Li. Then, Zj > x, if L > Li, and Zi < Xi \i L < L^. Also Zj = 1. Thus, 

~ 2^ Li ^ Li ^ L, - ^ L ^ L ■ 

i i:L>Li i:L<Li i:L>Li i:L<Li 

Moreover, V = Q \i and only if Zi = Xi for all i if and only if x is a unit vector. 
Consider next the function ^^QiLi. Then, 

5:Q.L. = A^L. = AI) = = ^ = L; 



hence, Yli Qi^i is decreasing. 

The function A DjLj = L • ^ ■ XjLj is the product of decreasing functions and hence 
decreasing. 

Finally, let V = Ei>2(^i In A - Li InDi). Then 

T/ StIt^' t sr(rQi-Di , Qi-Di \ 

^ E (l.^^^^^ - L,^i^Z^) ^ E(^^ - < 0- 



A Lli , 

The Lyapunov function ^•>2(Ljln A — LilnL^i) was already considered in |MU07j . 

Theorem 1 (Miyashi-Ohnishi [ MO07| ) For a network of parallel links, the dynamics 
converge against Di = 1 and A = for i > 2. 

Proof: xi = Di/D is monotonically increasing and bounded by 1. Hence, it converges. 
Assume that the limit x^ is less than one. Clearly, x^ > 0. For xi < x^, we have 1/L = 
Yli^i/Li < xl/Li + {l — xl)/L2- Moreover, for large enough t, xi > xJ/2 and D <2 (Lemma 
El), and hence, x'l > e for some e > 0. Thus, < 1 is impossible. I 
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Some of the Lyapunov functions have natural interpretations: QiLi is the total cost 
of the flow; (Yli -^i-^i) / "l^i -^i is the total hardware cost normalized by the total diameter, 
where a link of length L and diameter D has cost DL; and A Y^- DiLi is the potential 
difference between source and sink multiplied by total hardware cost. These functions are 
readily generalized to general networks by interpreting the summations as summations over 
all edges of the network. Our computer simulations showed that none of these functions is a 
Lyapunov function for general networks. 

However, ^ • Di can also be interpreted as the minimum capacity of a source-sink cut in 
a network where Di is the capacity of edge i. With this interpretation, (YliiDiL.i) / Y^^Di 
becomes 

minsec Cs ' 

where C is the set of all sq-si cuts and Cs is the capacity of the cut C. Our computer 
simulations suggested that this function may serve as a Lyapunov function for general graphs. 
We will see below that a slight modification is actually a Lyapunov function. 

5 Electrical Networks and Simple Facts 

In this section, we establish some more notation, review basic properties of electrical networks, 
and prove some simple facts. 

Each node v of the graph G has a potential that is a function of time. A potential 
difference Ae between the endpoints of an edge e induces a flow on the edge. For e = {u,v), 

Qe = DeAe/Le = De{pu - Pv)/Le = {Pu " Pv) / Re (4) 

is the flow across e in the direction from u to v. If Qe < 0, the flow is in the reverse direction. 
The potentials are such that there is flow conservation in every vertex except for sq and si 
and such that the net flow from sq to si is one, that is, for every vertex u, we have 

^ Qu,v = b{u), (5) 

v:{u,v)£E 

where 6(so) = 1 = — 6(si) and b{u) = for all other vertices u. After fixing one potential 
to an arbitrary value, say = 0, the other potentials are readily determined by solving a 
linear system. This means that each Qe can be expressed as a function of R only. 

For the main convergence proof, we will use some fundamental principles from the theory 
of electrical networks (for a complete treatment, see for example [Bol98, Chapters II, IX]). 

Thomson's Principle. The flow Q is uniquely determined as a feasible flow that minimizes 
the total energy dissipation Y2eReQe^ with Re = Le/De- In other words, for any flow x 
satisfying ([S]), 

Y,ReQl<Y.^^'^e- (6) 

e e 

KirchhofF's Theorem. For a graph G = {N, E) and an oriented edge e = (n, v) G E, let 

• Sp be the set of all spanning trees of G, and let 

• Sp(u, v) be the set of all spanning trees T of G, for which the oriented edge (n, v) lies 
on the unique path from sq to si in T. 
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For a set of trees S, define T{S) = J2TeS IleeT ^e/Le- Then, the current through the edge e 
is 

_ r(Sp(n,^)) - r(Sp(7;,-u)) 

Gronwall's Lemma. Let a, /3 € M and let x be a continuous differentiable real function on 
[0, oo). If ax{t) < x{t) < /3x(t) for ah t > 0, then 

x(0) e"* < x{t) < x(0) e'^* for ah t > 0. 

Proof: 

d X xe^^ — I3xe^* x{t) x(0) 

A similar calculation establishes x{t) > x(0)e"*. I 

The next lemma gives some properties that are easily derived from ([1]), Q, and ([5]). 
Recall that C is the set of sq-si cuts and Cs = J2ee5{s)^<^- ^^so, let Lmin = niiUeLe, 
Lmax = maXgLg, n = \N\, and m = \E\. 

Lemma 4 The following hold for any edge e E and any cut S & C: 

(i) IQel < 1. 
(ii) EeeS{{so}) \Qe\ = 1- 

(in) De{t) > DeiO) exp(-t) for all t, 

(iv) De{t) < 1 + {De{0) - 1) exp(-t) for all t. 

(v) Re > -Lmin/2 for all sufficiently large t. 

(vi) Csif) > 1 + {Cs{^) — 1) exp(— t) for all t, with equality if S = {sq}- 

(vii) C|s(j} — > 1 as i —> oo. 

(via) Orient the edges according to the direction of the flow. For sufficiently large t, there is 
a directed source-sink path, in which all edges have diameter at least l/2m. 

(ix) |Ae| < 2nmLjnax for all sufficiently large t. 

(x) De/D(, G [— 1, 2nmLjnax/-^min] for all sufficiently large t. 

Proof: 

(i) Since Q is a flow, it can be decomposed into sq-si flow paths and cycles. If \Qe\ > 1, 
since &(so) = 1, there exists a positive cycle in this decomposition, a contradiction to the 
existence of potential values at the nodes. The claim is also an immediate consequence 
of ©. 

(ii) It follows from equations (H]) and ([5]) that = max^p„, so Qs(y,v > for all {sq, v} G E, 
and Y.ee5{{so}) l'3e| = Y.eeS{{so}) = ^• 
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(iii) From the evolution equation ([T]), > — Dg- The claim follows by Gronwall's Lemma. 

(iv) \Qe\ < 1 for any edge e, so D(, < 1 — D(, from ([1]), and the claim follows as before. 

(v) From ()iv|l . <2 for all sufficiently large t, so Re = L^/De > for the same t's. 

(vi) = Eee5(S) -^e = Eee5(5)(l<3e| - -Dg) > 1 - Cs, with equality if 5 = {sq}- 

(vii) Follows by noting that the inequality in (jvi]) becomes tight for the cut {sq}, due to 

(viii) From ()vi|) . eventually C5 > 1/2 for all S" € C, so there is an edge of diameter at least 
l/2m in every cut. Thus, there is a sq-si path in which every edge has diameter at least 
l/2m. 

(ix) Consider a source-sink path in which every edge has diameter at least l/2m. By the 
total potential drop — is at most 2nmLmax- 

(x) D(,/De = {\Qe\ — D(,)/D(, = [Ael/ie " Ij and the bound follows from (|ix|). | 



6 Convergence 

We will prove convergence for general graphs. Throughout this section, we will assume that 
t is large enough for all the claims of Lemma S] requiring a sufficiently large t to hold. 

6.1 Properties of Equilibrium Points. 

Recall that D e is an equilibrium point, when Df, = for all e £ E, which by ([T]) is 
equivalent to Dg = \Qe\ for all e € E. 

Lemma 5 At an equilibrium point, minggc C5 = C|sq} = 1. 
Proof: 

ee5{S) eeS({so}) m 



Lemma 6 The equilibria are precisely the flows of value 1, in which all source-sink paths 
have the same length. If no two source-sink paths have the same length, the equilibria are 
precisely the simple source-sink paths. 

Proof: Let Q be a flow of value 1, in which all source-sink paths have the same length. 
We orient the edges such that Qe ^ for all e and show that D = Q \s an equilibrium 
point. Let Ei be the set of edges carrying positive flow, and let V\ be the set of vertices lying 
on a source-sink path consisting of edges in Ei. For f G Vi, set its potential to the length 
of the paths from v to si in {Vi,Ei); observe that all such paths have the same length by 
assumption. Let Q' be the electrical flow induced by the potentials and edge diameters. For 
any edge e = (u, f) € -Ei, we have Q'^, = D^^elE^ = = Qe- Thus, Q' = Q. For any edge 
e ^ El, we have Qe = = -Dg- We conclude that D is an equilibrium point. 
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Let D be an equilibrium point and let Qe be the corresponding current along edge e, 
where we orient the edges so that Qe ^ for all e £ E. Whenever > 0, we have 

= QeLe/De = bccausc of the equilibrium condition. Since all directed sq-si paths 
span the same potential difference, all directed paths from sq to si in {e & E : > 0} have 
the same length. Moreover, by Lemma [SJ min^-Cs = 1. Thus, D \s a. flow of value 1. | 

Let 8* be the set of flows of value one in the network of shortest source-sink paths. If the 
shortest source-sink path is unique, £* is a singleton, namely the flow of value one along the 
shortest source-sink path. 

6.2 The Convergence Process 

Lemma 7 LetW = (C{^„} - l)^. Then, W = -2W < 0, with equality iff C^^^} = 1. 
Proof: Let Cq = Cj^g} for short. Then, since X^gg^d^o}) \Qe\ = 1, 

W = 2{Co - 1) Yl - = 2(^0 - - ^o) = -2(Co - if < 0. 

eeS{{so}) ■ 



The following functions play a crucial role. Let C = min^gc Cs, and 

1 



Vs = t;- LeDe for each G C, 



V = max Vs + W, and 

sec 

Lemma 8 Let S be a minimum capacity cut at time t. Then, Vs{t) < —h{t). 

Proof: Let X be the characteristic vector of 5(5), that is, X(. = 1 if e € 5{S) and 
otherwise. Observe that Cs = C since 5 is a minimum capacity cut. We have 



< ^ E - ^ E ^-^e - ^ E + ^ E 

e e e e 

= -h. 
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The only inequality follows from Le = ReDe and Xe\Qe\ > 1, which holds because at least 
one unit current must cross S. I 



Lemma 9 Let f(t) = maxs^c fsit); where each fs is continuous and differentiable. If f{t) 
exists, then there is S € C such that f{t) = fs{t) and f{t) = fsit). 

Proof: Since C is finite, there is at least one S € C such that for each fixed 6 > 0, f{t + e) = 
fs{t+£) for infinitely many e with \e\ < 6. By continuity of / and fs, this implies f{t) = fs{t)- 
Moreover, since 

j.^ max5/ fs'jt + e) - maxy /g/(t) 
£-s>o e 

exists and is equal to f{t), any sequence ei,e2, ■ ■ ■ converging to zero has the property that 
maxg/ fs'it + Ei) - maxs' fs'{t) 



Si 



f{t) for i — )• oo. 



Taking (eJi^i to be a sequence converging to zero such that f{t + e,) = fs{t + e,) for all i, 
we obtain 

Lemma 10 V exists almost everywhere. IfV{t) exists, then V{t) < —h{t) — 2W{t) < 0, and 
V{t) = if and only if D(,{t) = for all e. 

Proof: V is Lipschitz-continuous since it is the maximum of a finite set of continuously 
differentiable functions. Since V is Lipschitz-continuous, the set of f s where V{t) does not 
exist has zero Lebesgue measure (see for example |CLSW98l Ch. 3], |MN92t Ch. 3]). When 
V{t) exists, we have V{t) = W{t) + Vs{t) for some S of minimum capacity (Lemma[9l). Then, 
Vit) < -h{t) - 2W{t) by Lemmas [7] and [3 

The fact that > is clear. We now show that h > 0. To this end, let F represent a 
maximum sq-si flow in an auxiliary network, having the same structure as G, and where the 
capacity on edge e is set equal to Dg- In other words, F is an sq-si flow satisfying \Fe\ < 
for all e € -E and having maximum value. By the max-flow min-cut theorem, this maximum 
value is equal to C = minggc Cs- But then, 

e 

rA'^v^ 1 

RpD^ 1 — — 2 ReDl 





1 , 




c ^ 






< 












< 










1 


< 






C72 



/ \ e / e 

^ (E ^el^e) (^E ^el^e) " ^ E ^^^e 



= 0, 

where we used the following inequalities: 
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- the Cauchy-Schwarz inequality Zei^e^lQeDiRe^D,) < {Y^eReQlY'^T^eReDlfl^; 

- Thomson's Principle ([6]) applied to the unit-value flows Q and F/C; Q is a minimum 
energy flow of unit value, while F/C is a feasible flow of unit value; 

- the fact that \F(,\ < -De for all e G -E. 

Finally, one can have /i = if and only if all the above inequalities are equalities, which 
implies that \Qe\ = \Fe\/C = D^jC for ah e. And, = iff EeG5{{so}) = ^ = 
EeG5({so}) IQel- So, /i = = iff IQel = for ah e. I 

The next lemma is a necessary technicality. 
Lemma 11 The junction t i— )• h{t) is Lips chitz- continuous. 

Proof: Since is continuous and bounded (by ([I])), is Lipschitz-continuous. Thus, it 
is enough to show that Qg is Lipschitz-continuous for all e. 

First, we claim that D(.{t+e) < {l + 2Ke)D(. for all e < l/4:K, where K = 2nmLinax/-Z^min- 
For if not, take 

e = mf{6 < 1/4K : D{t + 6) > {1 + 2K6)D{t)}, 

then e > (since De{t) < KDe{t) by LemmaS]) and, by continuity, De{t+£) > {l+2Ke)De{t). 
There must be t' e [t,t + e] such that De{t') = 2KDe{t). On the other hand, 

De{t') < KDeit') < K{1 + 2Ke)De{t) 

< K(l + 2K/AK)De{t) < 2KDe{t), 

which is a contradiction. Thus, D(,{t + e) < (1 + 2K£)De for all e < l/AK. Similarly, 
De{t + e) > (1 - 2Ke)De. 

Consider now a spanning tree T of G. Let 7t = HeeT-^eZ-^e- Then ^T{i + e) < (1 + 
2Ke)"7T(t) < (l + 4ni^e)7T(t) for sufficiently smaU e. Similarly, 7T(t+e) > (l-4ni^e)7T(t). 

By Kirchhoff 's Theorem, 

^ _ SrgSp(n,i,) 7T - ETgSp(i.,») 7T 
2L.TGSP 7T 

and plugging the bounds for ^xit + e)/^T{f) shows that Qe{t + e) = Qe{t){l + 0(e)), where 
the constant implicit in the O(-) notation does not depend on t. Since \Qe\ < 1, we obtain 
that \Qe(t + £) — Qe{t)\ < 0(1) •£, that is, Qe is Lipschitz-continuous, and this in turn implies 
the Lipschitz-continuity of h. I 



Lemma 12 \De — \Qe\\ converges to zero for all e £ E. 

Proof: Consider again the function h. We claim /i ^ as t — )■ oo. If not, there is e > and 
an infinite unbounded sequence ti,t2, ■ ■ ■ such that h{ti) > e for all i. Since h is Lipschitz- 
continuous (Lemma [TT|) . there is 6 such that h{ti + 6') > h{ti) — e/2 > e/2 for all 5' G [0,5] 
and all i. So by Lemma [TOl V{t) < —h{t) < —e/2 for every t in [ti^ti + 5\ (except possibly a 
zero measure set), meaning that V decreases by at least £6/2 infinitely many times. But this 
is impossible since V is positive and non-increasing. 
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Thus, for any e > 0, there is to such that h{t) < e for ah t > to- Then, recalhng that 
Re ^ -^min/2 for ah Sufficiently large t (Lemma IHvj) . we find 

e ^ ^ e ^ ^ 

e e e 

e e 

= 2/i < 2e, 

where we used once more the inequality "^^ReQe ^ J2e ^^D'^/C'^ , which was proved in 
Lemma 1101 This implies that for each e, D^/C — \Qe\ — ?> as t — )• oo. Summing across 
e G '^({•5o}) and using Lemma r4lii[ we obtain C^g^^/C — 1 ^ as t ^ oo. From Lemma [H 
C{sq} — )• 1 as t — )• oo, so C ^ 1 as well. 

To conclude, we show that D^/C — \Qe\ and C — > 1 together imply De — \Qe\ — > 0. 
Let e > be arbitrary. For all sufficiently large t, \De/C — \Qe\\ < e, |1 — C*] < e, < 2, 
and C > 1/2. Thus, 

\De - \Qe\\ < \De " D,/C\ + \DJC - \Qe\\ < dJ-^^ + \De/C - \Qe\\ < bs. _ 



Lemma 13 Let A = PsQ—Psi be the potential difference between source and sink. A converges 
to the length L* of a shortest source-sink path. 

Proof: Let £. be the set of lengths of simple source-sink paths. We first show that A 
converges to a point in C and then show convergence to L*. 

Orient edges according to the direction of the flow. By Lemma I4lviiil there is a directed 
source-sink path P of edges of diameter at least 1 /2m. Let e > be arbitrary. We will show 
|A — Lp\ < e. For this, it suffices to show |Ae — Le\ < e/n for any edge e of P, where Ag 
is the potential drop on e. By Ohm's law, the potential drop on e is Ag = {Qe/De)Le, and 
hence, |Ae — Le| = \Qe/De — l\Le = \{Qe — De)/ De\Le < 2mLmax | Qe " -Ce | • The claim follows 
since \Qe — De\ converges to zero. 

The set C is finite. Let e be positive and smaller than half the minimal distance between 
two elements in C By the preceeding paragraph, there is for all sufficiently large t a path Pt 
such that |A — I < e. Since A is a continuous function of time, Lp^ must become constant. 
We have now shown that A converges to an element in C. 

We will next show that A converges to L*. Assume otherwise, and let P' be a shortest 
undirected source-sink path. Let Wp' = YleeP' -^elu-De- This function was already used by 
Miyaji and Ohnishi [MU08] . We have 

^P'= E ^i\Qe\ - De) = J2 \^e\-J2^-^P^0-Psi-Lp' = A-L*. 
e&P' ef^pi eeP' 

Let 5 > be such that there is no source-sink path with length in the open interval {L*,L* + 
25). Then, A — L* > 6 for all sufficently large t, and hence, Wp' > 5 for all sufficiently large 
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t. Thus, Wp' goes to +00. However, Wp' < riL^s^x for all sufficiently large t since <2 for 
all e and t large enough. This is a contradiction. Thus, A converges to L*. | 

Lemma 14 Let e be any edge that does not lie on a shortest source-sink path. Then, Dg and 
Qe converge to zero. 

Proof: Since \D(, — \Qe\\ converges to zero, it suffices to prove that converges to zero. 
Assume otherwise. Then, there is a (5 > such that \Qe\ > 5 for arbitrarily large t. 

Consider any such t and orient the edges according to the direction of the flow at time 
t. Let e = («, f). Because of flow conservation, there must be an edge into u and an edge 
out of V carrying flow at least Qe/n. Continuing in this way, we obtain a source-sink path 
P in which every edge carries flow at least Qe/n^ > P depends on time and Lp > L* 

always. We will show |A — Lp\ < (Lp — L*)/4 for sufficiently large t, a contradiction to the 
fact that A converges to L*. For this, it suffices to show \Ag — Lg\ < {Lp — L*)/{An) for 
any edge g of P, where A^ is the potential drop on g. By Ohm's law, the potential drop 
on fit is Ag = {Qg/Dg)Lg, and hence, |Ag - Lg\ = \Qg/Dg - l\Lg = \{Qg - Dg)/Dg\Lg < 
Lra..\Qg " Dg\/Dg. For large enough t, \Qg - Dg\ < min(V(2n"), 5(Lp - L*)/(8n"+iL^ax)). 
Then, Dg>Qg- \Qg - Dg\ > V(2n"), and hence, L^,,\Qg - Dg\/Dg < {Lp - L*)/(4n). | 

Theorem 2 The dynamics are attracted by £* . Lf the shortest source-sink path is unique, 
the dynamics converge against a flow of value 1 on the shortest source sink path. 

Proof: Q is a source-sink flow of value one at all times. We show first that Q is attracted 
to E* . Orient the edges in the direction of the flow. We can decompose Q into fiowpaths. For 
an oriented path P, let lp be the unit flow along P. We can write Q = Y^pXplp, where xp 
is the flow along the path P. This decomposition is not unique. We group the flowpath into 
two sets, the paths running inside Gq and the paths using an edge outside Go, i.e., 

Q = <3o + Qi, where Qo = ^ xplp. 

P is a path in Go 

Qo is a flow in Go, and each flowpath in Qi is a non-shortest source-sink pathjl We show 
that the value of Qo converges to one. 

Assume otherwise. Then, there is a (5 > such that the value of Qi is at least 5 for 
arbitrarily large times t. At any such time, there is an edge e ^ Eq carrying flow at least 6/m; 
this holds since source-sink cuts contain at most m edges. Since there are only finitely many 
edges, there must be an edge e ^ Eq for which Qg does not converge to zero, a contradiction 
to Lemma UM 

We have now shown that the distance between Q and £* converges to zero. By Lemma [T2l 
\De — \Qe\\ converges to zero for all e, and hence, the distance between Q and D converges 
to zero. Thus, D is attracted by £*. 

^The decomposition into Qo and Qi can be constructed as follows: Initialize Qo to Q and Qi to the empty 
flow. Consider any edge e ^ Eq carrying positive flow in Qo, say e. Let P be an oriented source-sink path 
carrying e units of flow and using e. Add elp to Qi and subtract it from Qo. Continue until Qo is a flow in 
Go- 
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Finally, if the shortest source-sink path is unique, £* is a singleton, and hence, D converges 
to the flow of value one along the shortest source-sink path. I 



Lemma 15 // the shortest source-sink path is unique, converges to dist(f,si) for each 
node V on the shortest source-sink path, where dist(7;,si) is the shortest path distance from v 
to Sl. 

Proof: Let Pq be the shortest source-sink path. For any e G P, converges to one and 
\De — Qe\ converges to zero. Thus, converges to L^. I 



6.3 More on the Lyapunov Function V 

In this section, we study V = LgDe/C -|- {C^g^^y — 1)^ as a function of D. Recall that 
C = C{D) = minsec Cs, where Cs = J2e(^S{S) ^e- 

Lemma 16 Let and be two equilibrium points. Define 

= {I - X)D^ + XD\ AG [0,1]. 

IfV{D^) < V{D^), then V{D'^) is a linear, increasing function of X. 

Proof: By Lemma[5l C{D^) = C{D^) = 1, and = C{s^y{D^) = 1. Since Cs{D) 

is linear in D for any fixed cut S, one has Cs{D^) > 1 and Cs{D^) > 1, so Cs{D^) > 1 
for all S. Thus, C{D^) > 1. On the other hand, q,^|(Z)^) = 1. Thus, C{D^) = 1, and 
V{D^) = LeD^, that is, V{D^) is a linear function of D^. I 



Lemma 17 The problem of minimizing V{D) for D € is equivalent to the shortest path 
problem. 

Proof: By introducing an additional variable C = ming Cs > 0, the problem of minimizing 
V{D) is equivalently formulated as 

2 



\ee5({so}) / 



mm 



s.t. Cs>C eC 
C > 
D>0. 



Substituting Xe = D^/C, we obtain 



mm 

S.t. J]] Xe > 1 VS" G C 

ee5(S) 

X > 0,C > 0, 



^ 2 
1 ' 

c 
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which is easily seen to be equivalent to the (fractional) shortest path problem. 



Lemma [T71 was the basis for the generalization of our results to the transportation problem 
(Section E]). We first generalized the above Lemma to Lemma [551 and then used the Lyapunov 
function suggested by the generalization. 

7 Rate of Convergence for Stable Flow Directions 

The direction of the flow across an edge depends on the initial conditions and time. We 
do not know whether flow directions can change infinitely often or whether they become 
ultimately fixed. In this section, we assume that flow directions stabilize and explore the 
consequences of this assumption. We will be able to make quite precise statements about the 
convergence of the system. We assume uniqueness of the shortest source-sink path and add 
more non-degeneracy assumptions as we go along. 

An edge e = {u, v} h econiGs hoTizofitdl if lim^—^QQ \Pu - Pv\ = 0, and it becomes directed 
from u to V (directed from v to u) if pu > Pv for all large t (p^ > pu for all large t). An edge 
stabilizes if it either becomes horizontal or directed, and a network stabilizes if all its edges 
stabilize. If a network stabilizes, we partition its edges into a set Eh of horizonal edges and 
a set ^ of directed edges. If {u,v} becomes directed from u to v, then {u,v) G . 

We already know that the diameters of the edges on the shortest source-sink path (we 
assume uniqueness in this section) converge to one. The diameters of the edges outside Go 
converge to zero. The potential of a vertex v (z Gq converges to dist(v,si). For stabilizing 
networks, we can prove a lot more. In particular, we can predict the decay rates of edges, the 
limit potentials of the vertices, and for each edge the direction in which the flow will stabilize. 

Definition 1 (Decay Rate) Let r < 0. 

A quantity D{t) decays with rate at least r if for every e > there is a constant A > 
such that for all t 

D{t) < ^6^''+^)*, or equivalently, \nD{t) < {\nA) + (r + e)t. 

A quantity D{t) decays with rate at most r if for every e > there is a constant a > 
such that for all t 

D{t) > ae^^""^^*, or equivalently, In D{t) > (In a) + (r - s)t. 

A quantity D{t) decays with rate r if it decays with rate at least and at most r. 

We first establish a simple Lemma that, for any edge, connects the decay rate of the fiow 
across the edge and the diameter of the edge. 

Lemma 18 Let —1 < a < and let e,g E. If Qe decays with rate at least a, then so does 
De- De decays with rate at most —1. // ||Qe| — |Qg|| decays with rate at least a, then \De — Dg\ 
decays with rate at least a. 

Proof: Assume first that Qe decays with rate at least a, where —1 < a < 0. Then, for any 
£ > 0, there is an ^ > such that Qe < y4e('*+^)* for ah t. Consider / with / = Ae(''+^)* - /. 
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This has solution / = foe * + ae^""^^^^ , where a = A/{1 + a + e) and /o is determined by the 
value of / at zero, namely, /(O) = /o + a. Consider — f ■ Then, 

- /) = \Qe\ -D,- (^e('^+^)* - /) < -pe " /)• 

Thus, Df, — / < C'e^* for some constant C by Gronwall's Lemma, and hence, 

De < ifo + C')e-' + ae("+")* < C"e^''+'^' 

for some constant C" . Thus, decays with rate at least a. 

De = \Qe\ — Dg > —De- Thus, Dg decays with rate at most —1 by Gronwall's Lemma. 
Finally, assume that HQel — \Qf\ \ decays with rate at least a. Then, 

j^{De - Dg) = \Qe\ ~ \Qf\ - {De - Dg) < \\Qe\ " " (^e " Dg), 

and therefore, De — Dg decays with rate at least —a. The same argument applies to Dg — De- 



For a path P, let W{P) := Y^^^p LelnDe be its weighted sum of log diameters, and let 
A(P) = Pa — Pb be the potential difference between its endpoints. The function W{P) was 
introduced by Miyaji and Ohnishi |MO071 [MO08] . 

Lemma 19 Let P be an arbitrary path, let A(P) be the potential drop along P, and let 
W{P) = Y.eepLe\nDe. Then, 

W{P) = A{P)-L{P) + 2 

eeP: A(e)<0 

// A(P) < A and A(e) > —5 for some 5 > 0, all e G P and for all sufficiently large t, then 

W{P){t) < C + (A - L{P) + 2n6)t 
for some constant C and all t. If A{P) > A for all sufficiently large t, then 

W{P){t) >C+{A-L{P))t 

for some constant C and all t. 

Proof: The first claim follows immediately from the dynamics of the system. 

W{P) = |A(e)I - LiP) = A{P) - L{P) + 2 ^ |A(e)|. 

eeP e&P: A{e)<0 

Let to be such that A(P) < A and A(e) > —5 for all t > to. We integrate the equality 
from to to t and obtain 



W{P){t) - W{P){to) = [ W{P)dt < (A - L(P) + 2n6){t - to). 

Jtn 
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This establishes the claim for t > to- Choosing C sufficiently large extends the claim to all t. 
Let to be such that A(i-') > A. We integrate the equality from to to t and obtain 



W{P){t) - W{P){t,) 




W{P)dt > (A - L(P))(t - to). 



This establishes the claim for t > to- Choosing C sufficiently large extends the claim to all t. 

I 

Edges that do not lie on a source-sink path never carry any flow, and hence, their diameter 
evolves as Z)e(O) exp(— t). From now on, we may therefore assume that every edge of G lies 
on a source-sink path. 

Lemma 20 For e € E^, decays with rate —1, and \Qe\ decays with rate at least —1. 

Proof: We certainly have < 2 for all large t. Let e = {u,v}, and let e > be arbitrary. 
Then, \pu — Pv\ < ^-^e for all large t, and hence, \Qe\ = {De/ Le)\pu — Pv\ ^ ^De for all large 
t. Thus, De < {e — l)De for all large t, and hence, {d/dt)ln < —1 -|- £■ Thus, Dg decays 
with rate at least —1. Since D^. > -De, decays with rate at most —1. 

\Qe\ = {D(./ L(,)\pu —Pv\^ ADe for some constant A. Thus, \Qe\ decays with rate at least 
-1. I 

We define a decomposition of G into paths Pq to Pfc, an orientation of these paths, a slope 
f{Pi) for each Pi, a vertex labelling p*, and an edge labelling r. Po is 80 shortest sq-si path 
in G, f{Po) = 1, Te = /(-Pq) — 1 for all e G Po) and = dist(u,si) for all v G Pq, where 



where V is the set of all paths P in G with the following properties: 

- the startpoint a and the endpoint 6 of P lie on Pq U . . . U Pi_i, p* > pi, and /(P) = 



- no interior vertex of P lies on Pq U . . . U Pj-i; and 

- no edge of P belongs to Po U . . . U Pj-i. 

If p* > pI, we direct Pj from a to b. If p* = pi, we leave the edges in Pi undirected. We set 
Tg = /(Pj) — 1 for all edges of Pj, and p^ = pI + /(Pi)distp. {v, h) for every interior vertex v of 
Pj. Here, distp,(u,6) is the distance from v io h along path Pj. Figure [3] illustrates the path 
decomposition. 

Lemma 21 There is an i^ <k such that 



dist(7;,si) is the shortest path distance from v to si. For 




Pi = argmax/(P), 

Per 



{pl-pl)/L{P); 



f{Po) > /(Pi) > . . . > f{Pi,) > = /(Pio+i) 



f{Pk). 



''We assume that Po is unique. 

^We assume that Pi is unique except if f{Pi) — 0. 
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Figure 3: All edges are assumed to have length 1; Pq = (ei), Pi = (62,63,64), P2 = (65,66), 
P*so = 1' = 0, = 1/3, pl = 2/3, pl = 1/2, /(Pi) = 1/3, and /(P2) = 1/6. 
The path (62,65,65,64) has /-value 1/4. 

Proof: It suffices to show: if there is an i such that /(Pj+i) > f{Pi), then f{Pi) = /(Pj+i) = 
0. If no endpoint of Pj+i is an internal vertex of Pi, then /(Pj+i) = f{Pi); otherwise Pj+i 
would have been chosen instead of Pi. By assumption, equality is only possible if the /-values 
are zero. So we may assume that at least one endpoint of Pj+i is an internal node of Pi; call 
it c and assume w.l.o.g. that it is the startpoint of Pj+i. Split Pi at c into P^ and P?, and 
let d be the other endpoint of Pj+i; d may lay on Pi. 

Assume first that d does not lie on Pi and consider the path P^Pi^i. The /-value of this 
path is 

P*a-P*d ^ P*a-P*c+P*c- P*d 

LiPl) + L{Pi+i) L{Pl)+L{Pi^i)- 
Next, observe that (p* -pl)/L{Pl) = f{Pi) since p* is defined by linear interpolation and 
ip*c ~ Pd)/^(-f«+i) = f{Pi+i) ^ f{Pi)- III case of inequality, P/Pj+i is chosen instead of Pj. 
In case of equality, there are two paths with the same /-value. By assumption, this is only 
possible if the /-values are zero. 

Assume next that d also lies on Pj. We then split Pj into three paths P/, Pf, and Pf and 
consider the path P/Pj+iP?. We then argue as in the preceding paragraph. | 

Theorem 3 If a network stabilizes, then ^ = U,;<j„i?(Pj), the orientation of any edge 6 € ^ 
agrees with the orientation induced by the path decomposition, and Ef^ = Ui^if^E{Pi). The 
potential of each node v converges to p* . The diameter of each edge e € E \ Pq decays with 
rate r^. 

Proof: We use induction on i to prove: 

- for every vertex v £ Pq U . . . U Pi, the node potential p^ converges to ; 

- for every edge e € PqU . . .U Pmin(j,jo)' stabilizes in the direction of the path Pj 
containing e; 

- for every edge e G Pi L) ... U Pi, the diameter converges to zero with rate re, and the 
flow converges to zero with rate at leaslH r^. If 6 G Pj and i < io, tlie flow converges to 
zero with rate r^. 

^If for an edge e = {u,v}, Pu ~ Pv = always, then Qe = always. Thus, for horizontal edges, Qe may 
converge to zero faster than with rate —1. 
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Lemma [15] establishes the base of the induction, the case i = 0. Assume now that the 
induction hypothesis holds for i — 1; we establish it for i. Let P<i-i = PqU . . .U Pi-i- 
For e G E\P<i_i, let 

where Ve is the set of paths P' in G \ P<i-i from some a G P<i-i to some b G P<i~i 
with p* > and containing e. Then, maxg^p^. ^ /e = f{Pi)- For i < io, we have further 
f{Pi) > maxe^p<- fe > /(-Pj+i). In general, the last inequality may be strict; see Figure [3l 

Lemma 22 For e (z E\ P<j_i, \Qe\ and Df. decay with rate at least fe — 1. 

Proof: According to Lemma [THl it suffices to prove the decay of \Qe\- Let e £ E \ P<i-i 
and let e > be arbitrary. We need to show 

ln|Qe(t)| <C7 + (/e + e-l)t 

for some constant C and all sufficiently large t. 

If Qe{t) = 0, the inequality holds for any value of C. So assume Qe{t) ^ and also 
assume that the flow across e = {«, f} is in the direction from u to v. We construct a path 
R{t) containing uv. For every vertex, except for source and sink, we have flow conservation. 
Hence there is an edge {v,w) carrying a flow of at least Qe/n in the direction from v to w. 
Similarly, there is an edge (x, v) carrying a flow of at least Qe/n in the direction from x to v. 
Continuing in this way, we reach vertices in P<j_i. Any edge on the path R{t) carries a flow 
of at least Qe/n"'. 

Since potential differences are bounded by B := InniLj^gj^ (Lemma I4lixp . any edge e' on 
R{t) must have a diameter of at least QeLe/{n"'B) > {L^i^/ {n'"' B))Qe. Let c = Ljam/{n"B). 
Then, 

W{R{t))= Le'lnDe' > L{R{t)){lnc + ln\Qeit)\). 

The path R{t) depends on time. Let a{t) and h{t) be the endpoints of R{t). Since e does not 
belong to P<i-i, 

For large enough t, we have A{R{t)) < A*{R{t)) + eL{R)/2. Every edge e G R{t) either 
belongs to ^ or to Ey^ due to the assumption that the network stabilizes. In the former case, 
R must use e in the direction fixed in E , in the latter case, the potential difference across 
e converges to zero. We now invoke Lemma 1191 with 5 = £L[R) / {An). It guarantees the 
existence of a constant Ci such that 

W{R{t)){t) <Ci + {A*{R{t)) + eL{R)/2 - L{R) + eL{R)/2)t 

for all t. The constant Ci depends on the path R{t). Since there are only finitely many 
different paths we may use the same constant Ci for all paths R{t). 
Combining the estimates, we obtain, for all sufficiently large i, 

L{R{t)){\nc + In \Qe{t)\) < Ci + (A*(i?(t)) + eL{R{t)) - L{R{t)))t, 
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and hence, 

\n\Qe{t)\<Ci/L{R{t))-hic+{fe + e-l)t. 



Corollary 4 For e (z E \ P<j_i, \Qe\ and decay with rate at least f{Pi) — 1. If i < io, 
then for any e € -E \ P<i, \Qe\ and decay with rate at least f{Pi) — 5 — 1 for some (5 > 0. 

Proof: If i < iq, and hence, f{Pi) > 0, /e < f{Pi) for any edge e G E \ P<i. The claim 
follows. I 



Lemma 23 Let e € Pi. Then, D^, decays with rate f{Pi) — 1. If i < io, then \Qe\ decays with 
rate f{Pi) - 1. 

Proof: We distinguish the cases /(Pi) = and f{Pi) > 0. If /(Pi) = 0, the diameter of all 
edges e Pi decays with rate at least —1 (Lemma llOp . No diameter decays with a rate faster 
than —1. 

We turn to the case / := /(Pj) > 0. The flows across the edges in \ P<i decay with rate 
at least / — 1, and the flows across the edges edges in \ P<i decay faster, say with rate at 
least / — (5 — 1 for some positive S (Corollary [3D . We first show 

WiPi) <C + L{Pi) ■ max(lnZ)e, {f - 5 - l)t) (8) 

for sufficiently large t and some constant C. If Pi consists of a single edge e, W[Pi) = 
Leln De{t) and ([8]) holds. Assume next that Pi = ei . . . with k > 1. Consider any interior 
node u of the path. The flow into u is equal to the flow out of u, and u has two incident 
edge^lll in Pi. The flow on the other edges incident to u decays with rate at least f — 5 — 1. 
Thus for any two consecutive edges on Pj, | {Q^^ \ — \Qej+i \ \ decays with rate at least f — 6 — 1. 
By Lemma [T8| this implies that {De^ — -Dej+il decays with rate at least f — 6 — 1. Thus, we 
have = + gej, where \gej \ < Cie^'f~^~^^^ for some constant Ci and all j. Plugging into 
the definition of W{Pi) yields 

W{Pi)< ^ LeJn(2max(De,9e,)) 

< L(Pi) In 2 + L{Pi) max(lnP>e, In Cie^^"*^"^)*), 

and we have established dH]). 

Let to be large enough such that |A(Pi) - A*(Pi)| < 5L{P,)/2 for ah t > to. Then, by 
Lemma \T9\ 

W{Pi)>A + L{Pi){f-S/2-l)t (9) 

for some constant A and all t. 
Combining ([8]) and ([9]) yields 

A + L{Pi){f - 5/2 -l)t<C + L{Pi) ■ max(lnP)e, {f - 5 - l)t). 
'^Here, we need uniqueness of Pi. Otherwise we would have a network of paths with the same slope. 
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Thus, for every t we have either 

A + L{Pi){f - 5/2 -l)t<C + L{Pi) ■ InDe 



or 

A + L{P.i){f - 5/2 -l)t<C + L{Pi) -if -5- l)t. 

The latter inequahty does not hold for any sufficiently large t. Thus, the former inequality 
holds for all sufficiently large t, and hence, decays with rate at most f{Pi) — 1. By 
Lemma dHl \Qe\ cannot decay at a faster rate if /(Pj) > 0. I 



Lemma 24 For v Pi, the potentials converge to p*. For e E Pi and i < iq, the flow 
direction stabilizes in the direction of Pi . 

Proof: Assume i < io first. Let Pi = ei . . . e^- The flows and the diameters of the edges in 
Pi decay with rate f{Pi) — 1 (Lemma [23]). The flows and diameters of the edges incident to 
the interior vertices of and not on Pj decay faster, say with rate at least f (Pi) — 5 — 1, where 
5 > 0. For large t and any interior vertex of Pi, one edge of Pi must, therefore, carry flow into 
the vertex, and the other edge incident to the vertex must carry it out of the vertex. Thus, 
the edges in Pj must either all be directed in the direction of Pj or in the opposite direction. 
As current flows from higher to lower potential, they must be directed in the direction of Pj. 

Because the flow and the diameters of the edges not on Pj and incident to interior vertices 
decay faster, we have for any e > and sufficiently large t 

Qej = Qei (1 + £j) and P>e, = -Dei + 4)' 

where \ej\, \e'^ \ < e. The potential drop Ag^. on edge ej is equal to 
and hence, the potential drop along the path is 



p,-p, = Y^A,^ = ^m){l + e"), 



where e" goes to zero with e. The potential drop along the path converges to ~ Pi- 
Thus, Qe^/De^ converges to f{Pi), and therefore, the potential of any interior vertex v of Pj 
converges to p* . 

We turn to the case i > i^. The potentials of the endpoints of Pj converge to the same 
value. Thus, the potentials of all interior vertices of Pj converge to the common potential of 
the endpoints. I 

We have now completed the induction step. I 
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Figure 4: The Wheatstone graph. 

8 The Wheatstone Graph 

Do edge directions stabihze? We do not know. We know one graph class for which edge 
directions are unique, namely series-parallel graphs. The simplest graph which is not series- 
parallel is the Wheatstone graph shown in Figured! We use the following notation: We have 
edges a to e as shown in the figure. For an edge x, Rx = L^/Dx denotes its resistance and 
Cx = Dx/Lx denotes its conductance!! For edges a, 6, c, and d, the direction of the flow is 
always downwards. For the edge e, the direction of the flow depends on the conductances. 
We have an example where the direction of the flow across e changes twice. 

A shortest path from source to sink may have two essentially different shapes. It either 
uses e, or it does not. If e lies on a shortest path, Lemma [1^1 suffices to prove convergence as 
observed by [MOnSj . If (a, e, d) is a shortest pathl, let P = (a, e) and P' = (b). Then, 

-{W{P) - W{P')) > A(P) - L{P) - (A(P') - L{P')) = L{P') - L{P) > 0. 
dt 

Since W{P) is bounded, this implies W{P') — )■ — oo. Thus, Di, converges to zero. Similarly, 
Dd must converge to zero. More precisely, W{P') goes to — oo linearly, and hence, Df^ and 
similarly decay exponentially. 

The non-trivial case is that the shortest path does not use e. We may assume w.l.o.g. that 
the shortest path uses the edges a and c. The ratio 

Ra + Rc 1 + Rc/ Ra ^ + Ca/Cc Ca + Cc 

is the ratio of the resistance of a to the total resistance of the right path; define x^, Xc, and 
Xd analogously. Observe Xa + Xc = I and Xf, + = 1. Let 

X* — . 

La + Lc 

define x^, x*, and x^ analogously. Without edge e, the potential drop on the edge a is x^ 
times the potential difference between source and sink. If Da = Dc, which we expect in the 
limit, Xa = x*- 



Observe that we use the letter C with a different meaning than in preceding sections. 
^For simplicity, we assume uniqueness of the shortest path in this section. 
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Lemma 25 Let S = CaCbiCc + Q) + {Ca + Cb)CcCd + + Cb){Cc + Cd)Ce. Then, 

^"^^ ^{Cb + Cd + Ce){La + L,){Ca + C,){xl - Xa) + CeChL 



^H^d\yb + ^d) \ Xc 
Proof: The derivatives of Ca to Cg were computed by Miyaji and Ohnishi [MO07j : 

C 

Ca = -^-f—{ChCc + CcCd + CcCe + CdCe) — Ca 
Cc 

Cc = {CaCd + CaCii + CaCe + C^Ce) — Cc- 

The derivatives of Cb and Cd can be obtained from the above by symmetry (exchange a with 
b and c with d). We now compute Xa- 

d Cc —{CaCc — CaCc) 



dt Ca + Cc {Ca + Ccy 

- {-WtSCbCc + CcCd + CcCe + CdCc) - Ca) C, 



+ 



{Ca + CcY 

'stJ' 



c 

f 



Ca ( '§f~{CaCd + CaCi) + CaCg + C^Ce) — Cc 



{Ca + Cc)^ 

CaCc ( CaCd + CaCh + CaCe + CfjCg C^Cc + CcCd + CcCc + CdC^ 



S{Ca + Cc)^ \ -Lc -C'a 



SLaLc{Ca + Cc)^ 



((a + Crf + Ce)(Da - I^c) + Cc{CbLa - CdLc)) 



^^^^ ^ {Cb + Cd + Cc){Da - Dc) + CcCbLc 



SLaLc{Ca + CcY " ' - ' - ' ^ " '\Lc Ld/Dd 

Finahy, observe 

La Cc La{Ca + Cc) — Cc{La + Lc) Da — Dc 



Xa Xa 



La + Lc Ca + Cc {La+Lc){Ca + Cc) {La + Lc){Ca + Cc)' 



We draw the fohowing conclusions: 

- if Ce = 0, then sign(a;a) = sign(L'a — Dc) = sign(x* — Xa)- Thus, Xa converges mono- 
tonicahy against x*. 

- From + Xrf = 1 and x* + x* = 1, we conclude 

sign ( ^ - — ) = sign(x* - x;,). 
Vx* XdJ 
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Figure 5: The transition diagram under the assumption x* < x^. 

- if s = sign(a;* — Xb) = sign(x* — Xa), then sign(xa) = s. 

- if Xa,Xb > X*, then Xa decreases. 

- if Xa,Xb < X*, then Xa increases. 

- if X(i,Xc > x'^, then x^ decreases (equivalent to: if Xa,Xb < xl, then x^ increases). 

- if X(i,Xc < x'^, then Xd increases (equivalent to: if Xa,Xb > x^, then Xb decreases). 
Theorem 5 Assume < x^, that is, La/Lc < Lf^/Ld- Then, 

1. The regime Xa,Xb > x^ cannot be entered. By symmetry, the regime Xa,Xb < x* cannot 
be entered. 

2. In the regime Xa,Xb € Xa decreases and Xb increases. Hence, in this regime, the 
direction of the middle edge e can change at most once. 

3. If the dynamics stay in the regime Xa,Xb > xl forever, Xa and Xb converge. 

4. If the dynamics stay in the regime Xa,Xb < x* forever, Xa and Xb converge. 

Proof: At (1): In the regime Xa,Xb > x^, Xa and Xb both decrease, and hence, the dynamics 
cannot enter the regime from the outside. More precisely, we consider two cases: Xb > xl and 
Xa = xl, or Xa > and Xb = xl. 

If Xb > xl and Xa = xl, Xa is non-increasing, and hence, we cannot enter the regime. 

If Xa > xl and Xb = xl, Xb is non-increasing, and hence, we cannot enter the regime. 

At (2): Obvious from the equations. 

At (3): Then, Xa and Xb are monotonically decreasing and hence converging. The deriva- 
tive of Xb clearly goes to zero if x^ and Xa converge to x^. 

At (4): Symmetrically to (3). I 
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In Figure O we use S, M, and L to denote the three ranges: S = [0,x*], M = [x*,xl], 
and L = [xl,l]. The box M x M is divided into the triangles Xa < x^ and Xa > Xf,. The figure 
also shows that the boxes S x S and L x L cannot be entered and that the latter triangle 
cannot be entered from the former. 

We conclude the following dynamics: Either the process stays in S" x S or L x L forever 
or it does not do so. If it leaves these sets of states, it cannot return. Moreover, there is no 
transition from the set of states RL to the set of states LR. Thus, if the process does not stay 
in S* X S or L X L forever, the direction of the middle edge stabilizes. 

Assume now that the dynamics stay forever in S x S, or in L x L. Then, Xa and Xb 
converge. Let x^ and be the limit values. If the limit values are distinct, the direction of 
the middle edge stabilizes. If the limit values are the same, the edge is horizontal and hence 
stabilizes. We summarize the discussion. 

Theorem 6 The dynamics of the Wheatstone graph stabilize. 

9 The Uncapacitated Transportation Problem 

The uncapacitated transportation problem generalizes the shortest path problem. With each 
vertex v, a supply/demand hy is associated. It is assumed that Ylv^^ = 0. Nodes with 
positive are called supply nodes and nodes with negative 6^ are called demand nodes. In 
the shortest path problem, exactly two vertices have non-zero supply /demand. A feasible 
solution to the transportation problem is a flow F satisfying the mass balance constraints, 
i.e., for every vertex f , is equal to the net flow out of v. The cost of a solution is F^Le- 
The Physarum solver for the transportation problem is as follows: At any fixed time, the 
current Q is a feasible solution to the transportation problem satisfying Ohm's law (jH). The 
dynamics evolve according to ([1]). 

For technical reasons, we extend G by a vertex sq with 6^0 = connect sq to an arbitrary 
vertex and decrease b^ by one. The fiow on the edge (so,v) is equal to one at all times. 

Our convergence proof for the shortest path problem extends to the transportation prob- 
lem. A cut 5 is a set of vertices. The edge set 5{S) of the cut is the set of edges having 
exactly one endpoint in 5", and the capacity Cs of the cut is the sum of the D-values in the 
cut. The demand/supply of the cut is bs = YlveS^'"- ^ non-trivial if bs / 0. We 

use C to denote the family of non-trivial cuts. For a non-trivial cut S, let Fs = Cs/bs, and 
let F = min {Fs ; S* G C}. One may view F as a scale factor; our transportation problem has 
a solution in a network with edge capacities D(,/F. A cut S with Fs = F is called a most 
constraining cut. 

Properties of Equilibrium Points. Recall that D £ Rf is an equilibrium point when 
-De = for all e (z E, which is equivalent to = \Qe\ for all e (z E. 

Lemma 26 At an equilibrium point, mins^c Cs/\bs\ = C{so}/t>{so} = 1- 



Proof: 
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Lemma 27 The equilibria are precisely the solutions to the transportation problem with the 
following equal-length property: Orient the edges such that Qe >0 for all e, and let N be the 
subnetwork of edges carrying positive flow. Then, for any two vertices u and v, all directed 
paths from u to v have the same length. 

Proof: Let Q be a solution to the transportation problem satisying the equal-length prop- 
erty. We show that D = Q is an equilibrium point. In any connected component of N, fix 
the potential of an arbitrary vertex to zero and then extend the potential function to the 
other vertices by the rule Ag = Lg. By the equal- length property, the potential function is 
well defined. Let Q' be the electrical flow induced by the potentials and edge diameters. For 
any edge e = {u,v) £ N, we have Q'^ = DeA-e/Le = De = Qe- For any edge e ^ N, we have 
Qe = = De. Thus, D is an equilibrium point. 

Let D be an equilibrium point and let Qe be the corresponding current along edge e. 
Whenever De > 0, we have Ag = QeLe/De = Le because of the equilibrium condition. Since 
all directed paths between any two vertices span the same potential difference, A'^ satisfies the 
equal-length property. Moreover, by Lemma [26l ToinsCs/bs = 1, and hence, Q is a solution 
to the transportation problem with the equal-length property. I 

Let 8 be the set of equilibria and let 8* be the set of equilibria of minimum cost. 



Lemma 28 Let W = {C{so} - 1)^- Then, W = -2W < with equality iff C {so} = 1- 
Proof: Let Cq = C{s„} for short. Then, since Y.ee5({so}) \Qe\ = 1, 



W = 2(Co - 1) Yl (l^-l - = 2(^0 - 1)(1 - Co) = -2(Co - 1)2 < 0. 



eeS{{so}) 



The following functions play a crucial role: Let F = min^gc Fs, and 




V 



max Vs + W, and 

Sec 



h 




Lemma 29 Let S be a most constraining cut at time t. Then, Vs{t) < —h{t). 
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Proof: Let X be the characteristic vector of 5(5), that is, X(, = 1 if e G 5{S), and Xf. = 
otherwise. Observe that Fg = F since 5 is a most constraining cut. Let C = Cs- We have 

e 

= [leC -Y.L,,D,,xyi {\Qe\-D,) 

< ^ E ^elali^e - § E i?el?e^ - ^ ^ ^^^e + ^ ^ ^^^^ 

e e e e 

The only inequahty follows from Lg = ReF>e and ^^XelQel ^ l^sl) which holds because at 
least hs units of current must cross S. I 



Lemma 30 V exists almost everywhere. IfV{t) exists, then V{t) < —h{t) — 2W{t) < 0, and 
V{t) = Oiffye, De{t) = 0. 

Proof: The almost everywhere existence of V is shown as in Lemma [TOl 

The fact that > is clear. We now show that h > 0. To this end, let / represent a 
solution to the (capacitated) transportation problem in an auxiliary network having the same 
structure as G and where the capacity of edge e is set equal to Dg/F; f exists by Hoffman's 
circulation theorem |Sch031 Corollary 11. 2g]: observe that for any cut T, Ft > F, and hence, 
|6t| < Ct/F. Then, 

= I E Re\Qe\De - j^Yl ^^^e 

e e 

^ I (e ^^Q'^ (e ^^^e) - E ^^^e 

^ I (E ^^fe) (E ^el^e) - ^^^e 
^ ^ (E ^^^e) (E ^^^e) - ^ E ^^^e 

= 0, 

where we used the following inequalities: 

- the Cauchy-Schwarz inequality Ze(.Re^\Qe\)(.Re^D,) < (Ee ReQ^V^HEe ReD^)^/^; 
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- Thomson's Principle ([6]) applied to the flows Q and /; Q is a minimum energy flow 
solving the transportation problem, while / is a feasible solution; and 

- the fact that \fe\< De/F for all ee E. 

Finally, one can have /i = if and only if all the above inequalities are equalities, which 
implies that \Qe\ = \ fe\ = De/F for ah e. And, 1^ = iff Eee5({so}) = 1 = EeG5({so}) \Qe\- 
So, /i = = iff IQel = -De for ah e. I 



Lemma 31 The function t i— ?> h{t) is Lipschitz-continuous. 

Proof: The proof of Lemma [TT] carries over. I 
Lemma 32 \De — \Qe\\ converges to zero for all e ^ E. 

Proof: The first and last paragraph of the proof of Lemma [T2] carry over. We redo the 
second paragraph. 

The first paragraph establishes that for any e > 0, there is to such that h{t) < e for all 
t > tQ. Then, recalling that Re > Linin/2 for all sufficiently large t (by Lemma H]), we find 

e ^ ^ e ^ 

= E ^-^e + E ^-Qe - I E Re\Qe\De 

e e e 

<^,Y.^eDl-jY.^e\Qe\De 

e e 

= 2h< 2e, 

where we used once more the inequality Y^e^eQl < Ee^e^e/^^ which was proved in 
Lemma [30l This implies that for each e, De/F — \Qe\ — > as t — )• oo. Summing across 
e € 5({so}) and using Lemma |D|ii|) , we obtain C{s(,}/C — 1 — )• as t ^ oo. From LemmaHl 
C^so} — >■ 1 as t — 7> oo, so C ^ 1 as well. I 

We are now ready to prove that the set of equilibria is an attractor. 
Theorem 7 The dynamics are attracted by the set E of equilibria. 

Proof: Assume otherwise. Then, there is a network and initial conditions for which the 
dynamics has an accumulation point D that is not an equilibrium; such an accumulation 
point exists because the dynamics are eventually confined to a compact set. Let Q be the 
fiow corresponding to D. Since D is not an equilibrium, there is an edge e with De 7^ |Qe|- 
This contradicts the fact that \De — \Qe\ \ converges to zero for all e. I 



Theorem 8 // no two equilibria have the same cost, the dynamics converge to a minimum 
cost solution. 
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Proof: Consider any equilibrium D* , and let Q* = D* be the corresponding flow. Let T* 
be the edges carrying non-zero flow. T* must be a forest, as otherwise, there would be two 
equilibria with the same cost. Consider any edge e = {u, v) of T*, and let S be the connected 
component of T*\e containing u. Then Q* = b{S), and hence, distinct equilibria have distinct 
associated forests. We conclude that the set of equilibria is flnite. 

The y-value of D* is equal to the cost LeQl of the corresponding flow since W = 
and F = 1 in an equilibrium. If no two equilibria have the same cost, the F-values of distinct 
equilibria are distinct. 

y is a decreasing function and hence converges. Since the dynamics are attracted to the 
set of equilibria, V must converge to the cost of an equilibrium. Since the equilibria are a 
discrete set, the dynamics must converge to some equilibrium. Call it D*. 

We next show that D* is a minimum cost solution to the transportation problem. Orient 
the edges in the direction of the flow Q* . If Q* is not a minimum cost flow, there is an 
oriented path P from a supply node u to a demand node v such that Qe > for all edges of 
P, and P is not a shortest path from u to v. The potential difference A^^, converges to Lp. 
We now derive a contradiction as in the proof of Lemma [T3j 

Let P' be a shortest path from u to w in G, let L* = Lpi be its length, and let Wpi = 
^ggp/ Lg In De- We have 

Wp> = ^ ^(IQel - De) = \Ae\- Le>Pu- Pv - Lpi = Auv - L* . 

eeP' ^ eGP' eeP' 

Let 5 > be such that there is no path from n to f with length in the open interval (L*, L* + 
26). Then, A — L* > 6 for all sufficently large t, and hence, Wp' > 5 for all sufficiently large 
t. Thus, Wp' goes to +oo. However, Wp' < riL^s^x for all sufliciently large t since Df, <2 for 
all e and t large enough. This is a contradiction. I 



Lemma 33 The problem of minimizing V{D) for D € is equivalent to the transportation 
problem. 

Proof: By introducing an additional variable F = vums Cs /\b{S)\ > 0, the problem of 
minimizing V{D) is equivalently formulated as 

s.t. Cs/\b{s)\ >F ys eC 

F >0 
D>0. 

Substituting Xe = D^/F, we obtain 



mm 

e 



min LeXe + F^^'^ Xg 

\eG<5({so}) 

S.t. Y xe> \b{S)\ ysec 

x>0,F>0 



^ 2 
1 ' 

F 
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which is easily seen to be equivalent to the (fractional) transportation problem. 
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